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Dynamics of distributed sources is described by nonlinear partial differ- 
ential equations. Lagrangian analytical solutions of these (and associated) 
equations are obtained and discussed in the context of the Lagrangian mod- 
eling - from the Lagrangian Invariants to dynamics. Possible applications of 
distributed sources and sinks to geophysical fluid dynamics and to gravitation 
are indicated. 

Singular vortices, sources (sinks) and vortex-sinks are well known models 
in fluid dynamics (including geophysical fluid dynamics), magnetized plasma, 
superfluidity and superconductivity (see, for example, Refs. 1-3 and refer- 
ences therein). Dynamics of distributed vortices is well developed classical 
area of research. However, to our knowledge, dynamics of distributed sources 
has not been considered before. Structure and intensity of distributed sources 
can change in time and, in this sense, they are more rich and flexible models 
for various media than singular sources with fixed intensity. In this letter 
we show that distributed sources have a nontrivial dynamics, which in some 
cases can be described analytically in Lagrangian variables. 

Local intensity of distributed sources (sinks) in moving continuous media 
is characterized by the divergency of the velocity field a(t, x) = dvi{t, x.)/dxi, 
where summation over the repeated subscripts is assumed from 1 to the di- 
mension of the space s = 1,2,3. To get an equation for a , we consider 
superposition of localized sources, which move each other with induced ve- 
locity fieldt^''^^ Thus, the equation can be written in form: 

da da da . 

Here / may include diffusion and some other effects (see below). We start 
with consideration of "free" sources with / = 0. In this case a is the La- 
grangian Invariant (LI), which is conserved along trajectories. For the po- 
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tential motion with velocity Vi — d(p/dxi, from (1) we get: 

da _ dAif ^ dip dAcp ^ A A ~ (2) 
dt dt dxk dxk ' ' dx^ 

The moving media, described by equation (2), can be interpreted as the 
barotropic (or incompressible) fluid: dvi/dt — —dw/dxi with w determined 
by the equation Aw = —dvi/dxkdvk/dxi = —{d'^ip/dxidxk)'^. The equation 
for w is similar to the equation for the kinematic pressure in the incompress- 
ible fluid, but now the compressibility depends on the mechanism of the mass 
production in distributed sources. The evolution of the density p is deter- 
mined by the continuity equation: dp/dt — —ap -|- where p is the rate 
of mass production. One natural possibility is that source is self-adjusting 
in such a way that p = ap and fluid is incompressible. Let us stress that 
(2) is not equation for ordinary fluid, but a model equation for distributed 
sources. The physical nature of these sources can be different for different 
kind of media. Some possible applications of distributed sources and sinks 
are indicated below. Independently of an inducing mechanism of motion, if 
the source is distributed, we can expect that various parts of it will move 
each other. This simple interaction is reflected in equation (2). Dynamics 
of distributed sources, potentially, may depend on other effects, including 
deformation. Some generalizations of equation (2) are presented below. 

The nonlinear partial differential equation (2), generally, is not easy to 
solve, even numerically. The type of nonlinearity in (2) is the same as in the 
equation for two-dimensional (2D) vorticity fleld u = Aijj in incompressible 
fluid {ip is the stream function). However, there is a difference in the orienta- 
tion of the components of velocity, been expressed in terms of ^l) (instead of 
</?). We will see below that this leads to essential difference in dynamics and, 
in some sense, makes dynamics of 2D distributed sources more complicated 
than 2D vortex dynamics. It is known that 2D vortex dynamics, generally, 
is nonintegrable and chaotic [^'*^^. We can expect the same for 2D distributed 
sources in general situation. To obtain some analytical solutions of (2), we 
will use the fact that a is LI . 

For 2D axisymmetrical flow we have: 

, , 1 dirv) , , dr , , 

a{t,r)^ S — = "0(^0), -77 = (3) 

r or dt 

where v — dtp/dr is the radial velocity, subscript "o" indicate initial value and 
in the second equation in (3) r(t. To) is the trajectory. Introducing Lagrangian 
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variable p{t,Po) = r^(t, Tq), Po — r^, we rewrite (3) in the form: ^(^) = 
ao{ro) = P{Po), where we used that d/dt = d/dt . Multiphcation by 
dp/dpo, gives: g^(^) = f (^) = l^iPo)^- Solution of this equation has 
the form: dp/dpo = exp[P{po)t]. Integration of this solution with natural 
condition p = when po = 0,and time differentiation gives: 

p{t,Po) ^r\t,ro) ^ exp[P{z)t]dz,2rv{t,ro) = p{z) e^[P{z)t]dz 
Jo Jo 

(4) 

Formulas (4) present Lagrangian description of the system. The same 
formulas give Eulerian velocity v{t, r) and r) = ao{ro) in parametric form 
{to - parameter). We note that a differential rotation of the system will not 
affect the obtained results. The solution depends on the initial distribution 
ao{ro)- For example, let ao{ro) — 6 — const for Tq < I and ao{ro) = for 
r„ > /. From (4) in the Lagrangian description we get r{t,ro) = roexp{6t/2), 
v{t,ro) = {ero/2)exp{et/2) for < I and r{t,ro) = [l^ exp{9t) + - /2]V2^ 
v{t,ro) = l9l'^/2)exp{et)[fexp{9t) + - /2]-i/2 for > /. the Eulerian 
description v{t,r) = 9r/2 and a{t,r) — 9 ior r < L{t) — lexp{9t/2), for 
r > L{t) we have v{t, r) = {9f /2r) exp{9t) and a{t, r) = 0. 

The main thing is that distributed source {9 > 0) is spreading, while 
distributed sink {9 < 0) is shrinking. In contrast, for 2D incompressible 
fluid, where the role of LI plays the vorticity lv, axisymmetric flow is the 
stationary solution. Distributed 2D sinks and sources potentially can model 
some convective systems in the atmosphere. Particularly, such modeling can 
be useful for horizontally converging updrafts, leading to concentration of 
vorticity (compare with Ref. 2 and references therein), and for horizontally 
divergent downdrafts, which cause aircraft accidents (see Ref. 5 and refer- 
ences therein). 

For the tree-dimensional (3D) distributed sources with spherical symme- 
try, instead of (3), we have: a(t, r) = ^-§^.{r'^''--) — <^o(^"o); ^ = Introducing 
q{t,qo) = r^(t, To), Qo — and 7(^0) = cto(^o)) we get analogous solution: 



q(t,qo) =r (t,ro) = / exp[y(z)t]dz,3r v(t,ro) = / -f(z) exp[y(z)t]dz 
Jo Jo 

(5) 

For the one- dimensional (ID) distributed sources we will get general so- 
lution without any restriction of symmetry. First steps of obtaining solution 
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are the same as before for 2D and 3D: 

dv dx dx 

a{t,x) = ^ = (^o{xo), = '"^ = exp[ao{xo)t] (6) 

At this step we can not assume that x{t, Xg) — when Xo — 0, because, 
generally, we do not have symmetry relative to the origin of the coordinate. 
Instead, we have to impose a boundary condition in order to fix the velocity of 
the system of observation. One possibility is condition: v{t, oo) + v(t, — oo) = 
0. We will see that with symmetric initial distribution ao{—Xo) = ao{xo), 
this condition insures that x{t, 0) = 0. Time derivative of the last equation 
in (6) and integration over Xo gives: 

-Xo 



v(t,Xo) = / ao{z) exp[ao{z)t]dz + u(t),u(t) = v(t,0) (7) 
Jo 

Using indicated boundary condition, we have: 

'"(^) = / {ao{-z)exp[ao{-z)t] - ao{z)exp[ao{z)t]}dz (8) 
z Jo 

We see that u{t) = for symmetric initial distribution, as expected. Inte- 
grating (7) with respect to time, we get: 

rXo rt 

x{t,Xo) — Xo+ / {ex.p[ao{z)t\ — l}dz + / u{T)dT (9) 
Jo Jo 

Formulas (7) - (9) give general analytical description of the dynamics 
of distributed sources. With the simplest symmetric initial distribution 
Oi-oi^o) = ^ for I |< / and ao{xo) = for | |> we get solution 
analogous to the presented above for the axisymmetric 2D flow. Descrip- 
tion of interaction of several patches of distributed sources and sinks will be 
presented elsewhere. 

There is an analogy between distributed sinks and the gravitation, for 
which we also have general ID analytical solution in the Lagrangian variables^^l . 
Relative accelerations of gravitational masses depend on distances (in any di- 
mension) similarly to relative velocities of sinks^^'^]. Homogeneous solutions 
of the general relativity with the cosmological constant (CC)[^^ behave simi- 
larly to (5) with 7 = const. Incorporation of the 3D distributed sources into 
cosmological modeling can contribute to an explanation of the accelerated 
expansion of the Universe, which is now a hot topic in physics (see Ref. 8 
and references therein). 
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Let us consider a generalization of the equation for distributed sources: 
dJa_dJ„ dJ^ _dv da 

Here a = Xoit^x) is the inverse trajectory and a is an arbitrary parame- 
ter. With (7 = we return to (6). With cr = 1 we get ID analog of 3D 
vortex dynamics, for which Ukdai/dxk is the LI. This LI corresponds to 
the Kelvin's theorem of conservation of the velocity circulation in the lo- 
cal (Cauchy) form'^'^l It potentially can lead to formation of singularities 
due to the vortex stretching. With a — —1 equation (10) gives the inertial 
motion: v{t,a) — Vo{a). Prom (10) we get equation for the deformation: 
|(|£) = Q;o(a)(ff)^+". Solution of this equation is: |f = [1 - (Tta,{a)]-^. 
Here we assume that a ^ 0, the limit with a ^ gives (6). From (10) we 
now have: 

dv 

— = ao{a)[l - atao{a)]~^ (11) 

We see that for o" < 0, which includes inertial motion, we get finite-time 
singularity with compression («„ <0). For a > 0, which includes analog of 3D 
vortex dynamics, we get finite-time singularity with stretching {ag > 0). The 
intermediate case (7) - (9) has both tendencies (shrinking sinks and spreading 
sources), but avoids finite-time singularities. Wc note, that equation (10) can 
be written in the form (1) with / = cra'^ . Similar generalizations can be 
made for 2D and 3D distributed sources and can be interpreted in terms of 
the barotropic (or incompressible) fiuid with corresponding modification of 
w. 

Finally, let us try one more step in generalization by using second order 
spatial derivatives in LI: 

d ,d'^v^ d d d , , 

Assuming, for simplicity, boundary condition v{t, — oo) = and integrating 
(12) from — cxo to x, we get: ^ = |a^,Q; = Solution of this equation 
corresponds to (11) with a = 1/2. Surprisingly, we get equivalence of two 
LI: d'^v/dx^ in (12) and J1/2 = dv / dx{da/ dxY^'^ui (10). 

A more detailed analysis and interpretation of the obtained in this let- 
ter nonlinear partial differential equations, as well as their generalization 
with diffusion, will be presented elsewhere. Here we used them for illustra- 
tion of the Lagrangian modeling - from LI to dynamics - and corresponding 
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technique. This approach creates a broad variety of Lagrangian models for 
nonhnear continuous media. 
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